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Abstract 

The effective mass one-dimensional Schrodinger equation for the generaUzed Morse potential is 
solved by using Nikiforov-Uvarov method. Energy eigenvalues and corresponding eigenfunctions 
are computed analytically. The results are also reduced to the case of constant mass. Energy 
eigenvalues are computed numerically for some diatomic molecules. The results are in agreement 
with the ones obtained before. 
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I. INTRODUCTION 



In recent years, the effects of the coordinate-dependence of the mass on the solutions of the 
relativistic and/or non-relativistic wave equations have been received a great attentions [1- 
14]. This is so because it has a widely range apphcations on different areas, for example, 
in the study of the semiconductors [15], or of electronic properties of quantum wells and 
quantum dots [16], and the impurities in crystals [17-19]. Point canonical transformation 
[9], deformed algebras [20], quadratic algebra method [21], path integral method [22], Lie 
algebra approach [23-25], group-theoretical approach [26], and supersymmetric formalism 
[27] are some of the methods used in the literature to solve the wave equations for the case 
of constant and/or position-dependent mass (PDM) distributions. 

In the present work, we study the effects of the PDM on the solutions of the Schrodinger 
equation (SE) for the generalized Morse potential given by 

V{x) = Vie~^^^^'^^ -V2e~^^^'^\ (1) 

where Vi , and V2 are two parameters, and correspond to D, and 2D in the usual Morse po- 
tential, respectively, where D is the dissociation energy of the molecule, and the parameter 
/3 is the depth of the potential. The Morse potential is one of the important exponentially 
varying potential which describes the interaction in diatomic molecules [28] . The Morse po- 
tential is exactly solvable for s-waves. However, for any ^-states only numerical solutions can 
be obtained with different approximation techniques [29] . We intend to apply the Nikiforov- 
Uvarov (NU) method [30] to find out the energy spectra and corresponding eigenfunctions of 
the generalized Morse potential by using a suitable coordinate-dependence mass distribution 
function which is finite at infinity, and also enables us to solve the SE exactly. 

The work is organized as follows. In Section II, we find out the energy eigenvalues and 
corresponding eigenfunctions of the SE by using the NU-method in the case of PDM. We give 
also the results for the case of constant mass. We give the numerical results for the bound 
state energies of H2, and LiH molecules in the case of constant and spatially dependent 
mass in Table I. We summarize our conclusions in Section III. 
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II. NIKIFOROV-UVAROV METHOD AND BOUND STATES 



The one-dimensional time-independent effective mass SE equation with the potential 
V{x) reads [6] 



2m{x) dx^ 2 \m'^{x) I dx 



V{x) + Ua^^{x)-E^,j;{x)^0, 



(2) 



where prime denotes the derivative of mass m(x) with respect to position and 



{a + "f — a)m{x)m" (x) + 2(a — — a — ^)m''^{x) 



(3) 



Am^{x){a + 1) 

where the ambiguity parameters proposed by different authors in the literature satisfy a + 
f3 + J — —1. In the ordering proposed by Weyl, the parameters have the values a — 1, 
a — J — 0. In the parameter set due to Li and Kuhn a — a — 0, and 7 = — |, etc. [6]. It is 
found in Ref. [6] that the Weyl and Li and Kuhn ambiguity orderings are equivalent in the 
case of SE equation. 

By redefinition of the wave function 



ip{x) — \/ m{x)(j){x), 



we get the SE as 



(4) 



d^(f){x) 



2m{x) dx"^ 
The effective potential in the above equation is 



+ {Ueffix) - E)cf>{x) = 0. 



(5) 



U,ff{x) = Uo,p^{x) + V{x) + 



( 



3m'^(x 



— m 



(6) 



4m'^{x) \ 2m{x) 

We consider the generalized Morse potential, which can be used to describe the vibrations 
of a two- atomic molecule, as the following 



V{x) = Vie-'^^'' - Fse"^^ , (0 < x < oo) . 



(7) 
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where Vi, and V2 are two real parameters, and x — (r — ro). The parameter (3 is a' fro, here, 
a' is a positive parameter, and ro is the equihbrium distance [29]. 
Here, we prefer to use the following mass-distribution 



m(x) = 



UlQ 



(8) 



where mo corresponds to the constant mass, and the free parameter 77 satisfies the condition 
that 77 < 1 to handle a positive mass at the zero. We get a finite mass distribution with 
this constraint on rj. We can check out the results of the case of constant mass. Further, we 

solve the SE analytically by using such a mass-function. 
Substituting Eqs. (8) and (7) into Eq. (5) we get 



+ 



+ 



(1 - rje-l^'^y 



2moV2 



4/3V(A + A2) 



2moVi' 



(9) 



where 



A2 



a — a'y — a — 7 

1 + a 
q; + 7 — a 1 

2(1 + a) ^ 2 ■ 



3 

4' 



By using the transformation z = e {0 < z < 1), we obtain 



where 



dz'^ z{l — Tjz) dz [77(1 — rjz)]'' 



- Siz^ - e-iz - erii\4>{^) = , 



(10) 



- £1 
-£2 

—^ni 



V 



2q;7 — q; — 7 1 
2(1 + a) 4 
a + 7 + l\ , 2moV2 



1 + a 
2mo-B 



2?7lo^l 



(12) 
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Now to apply the NU-method [30], we rewrite Eq. (11) in the following form 



<t^"{z) + 44 + ^ m = 0, (13) 

where a{z) and a{z) are polynomials with second-degree, at most, and f{z) is a polynomial 
with first-degree. We define the total wave function as 



<t>{z) = azmz). (14) 
Thus Eq. (13) is reduced to a hypergeometric type equation 

a{z)tl)"{z) + t{z)^1]\z) + \tl){z) = 0. (15) 
We also define the new eigenvalue for the Eq. (13) as 



A = Xn^-nr'- (^^0,1,2,...) (16) 

where 



t{z) = f{z) + 2n{z). (17) 

The derivative of t{z) must be negative. A(A„) is obtained from a particular solution of 
the polynomial ipni^) with the degree of n. ipni^) is the hypergeometric type function whose 
solutions are given by [30] 



where the weight function p{z) satisfies the equation 



^^[a{z)p{z)]=T{z)p{z). (19) 
On the other hand, the function ^{z) satisfies the relation 



i\z)li{z)=7:{z)/a{z). 
Comparing Eq. (11) with Eq. (13), we have 



(20) 



T{z) = l-riz, a{z) = z{l - rjz) , a{z) = -Siz^ - e2Z - (21) 
The 7r(z) has the form [30] 

, , cr'iz) — fiz) , Lcr'(z) — f(z).^ ^, , , 7T 
^(^) = 2 ^ V 2 ^ " "^'^ + ^""^^^ ' ^^^^ 

or, exphcitly 



T^l-^) = - Y ± V ( J - + + (£2 + A;);^ + . (23) 

The constant is determined by imposing a condition such that the discriminant under 

the square root should be zero. The roots of k are ki^2 = —^2 — '^V^ni T \f^i-^-, where 
A = ^Ve„^ + 4r/e2 + r/2 + 4ei. 

Substituting these values into Eq. (23), we get for ki 



7t{z) = 



-qz 

T 



(I 



(24) 



and for k2 



A 



(25) 



Now we calculate the polynomial t(2;) from 77(2;) such that its derivative with respect to 
z must be negative. Thus we take the root k^-, and by using Eq. (25) for t^{z) we get 



(A \ 
t{z) = 1 + ly/e^i ~ ^ ( 2" "n^f^i + ?7 j ^ ■ 

with the derivative t'{z) — —2^^ + ri^Je^i + 77^ . 



(26) 



The constant X — k + n'{z) becomes 



A = -£2 - 2r)ene - V^t [A + r]) - ^{A + r]) , (27) 



and Eq. (16) gives us 



/'A \ 
= 2ni^- + 7jy/^e + Tjj + 7]n{n - 1) . (28) 

Substituting the values of the parameters given by Eq. (12), and setting A = A„, one can 
find the energy eigenvalues as 



^ ^ , [n' + n- {l/2)]rj - - 2[n + (1/2)] - (77V2) 

^ (2n + l)77- V4£i- 2772 ^ ' ^ ' 

It is seen that there is a strong dependence of the energy eigenvalues to the parameter rj. 
Further, one can see that the generalized Morse potential has a real energy spectra in the 
case of PDM under the condition that 



moV,^^,^a-2on-a-^_l_y ^^^^ 



f3^h' ' \ 1 + a 
We get the following energy spectra in constant mass case as 



2n + l + 



£2 



2 



(31) 



/ei- 

It is exactly same result obtained in the hterature [29] . 

We hst the numerical results for the bound state energies of the H2, and LiH molecules 
in the constant mass case in Table I. We use the same parameters given in Ref. (29), such 
as D, To, mo, a', and Eq, to compare our results. Here, the new parameter Eq is a short 
notation, i.e., fi^ /{rnQr"^). We also give the bound state energies of the above molecules for 
three different values of r] for each molecule in Table I. For simplicity, we choose the Weyl 
ordering for the ambiguity parameters in the numerical analyze. 

Now let us find the eigenf unctions. We first compute the weight function from Eqs. (17) 
and (19) 



p[z) = (1 _ r^zr , 



(32) 



where /I = ^/l + Aent + ^ (^2 + ^), and the wave functions become 



^n(^) 



[l-r,z) 



n+A 



(33) 



where &„ is a normahzation constant. In the hmit 77 — > 1, the polynomial solutions can be 
written in terms of the Jacobi polynomials with weight function p{z) as [31] 



{2r]z-l), A>-1, -2^i>-l. 



(34) 



On the other hand, the other part of the wave function is obtained from the Eq. (20) as 



^(^) =^-^(l-,y^)(V2)(l+A) 



Thus, the total eigenfunctions take 



" > 2 



'1 — X 



(35) 



</.„(^) = 6; Z-^' (1 - nz)^m{l+A)piA,-2^,) ^2^^ _ _ ^gg^ 

where h'^ is the new normalization constant. The normalization condition \(j)[z)\^dz = 1 
gives us 



.^^p04,-2V^) ^^^^^ ^ ^ ^ ^3^^ 



where we use a new variable defined as x = 2r]z — 1 . By using the required identities of the 
Jacobi polynomials [32, 33], we obtain 



X 



22A-4Vi;^+i 2{A- y/^ + 2n)+An{l + n) 
A - 2y/^ + 2n + l {A- 2^^ + 2n + 2)(i - 2^/^ + 2n) 
T{A + n+ l)r{-2^^ + n + l) 

n\r{A-2y/^ + n + l) 



(38) 
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III. CONCLUSION 



We have solved the one-dimensional position-dependent effective mass SE for the general- 
ized Morse potential by using NU-method and obtained analytically the energy eigenvalues 
and corresponding eigenf unctions. They depend on the free parameter r) strongly. We have 
shown that the results can be reduced to the ones obtained for the case of the constant 
mass. We have listed the numerical values of the energy eigenvalues in Table I for the H2, 
and LiH molecules for different values of the quantum number n, and free parameter 77. 
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TABLE I: The dependence of the bound states on n in eV for H2 {D = 4, 7446ey, ro = 0, 7416A, 
mo = 0,50391amM, a' = 1,440558, and Eq = 1,508343932 x IQ-^ey), and LiH molecules {D = 
2,515287eF, ro = 1,59561, mo = 0, 8801221am«, a' = 1,7998368, and Eq = 1,865528199 x 
10-=^eF) [29]. 



r/ = 



n 


En{H2r 








EniLiHY 


EniLiHY 





-4.476 


-4A76 






-2.429 


-2.429 


2 


-3.480 


-3.480 






-2.098 


-2.098 


4 


-2.609 


-2.609 






-1.792 


-1.792 


10 


-0.748 


-0.748 






-1.018 


-1.018 


15 


-0.057 








-0.539 




20 










-0.211 


-0.211 




n 


En{H2Y 


En{H2Y 


En{H2Y 


En{LiHY 


EniLiHY 


EniLiHY 





-4.528 


-4.582 


-4.637 


-2.446 


-2.463 


-2.481 


2 


-3.706 


-3.955 


-4.228 


-2.176 


-2.259 


-2.346 


4 


-2.953 


-3.363 


-3.856 


-1.920 


-2.062 


-2.219 


6 


-2.274 


-2.809 


-3.522 


-1.677 


-1.872 


-2.099 


10 


-1.152 


-1.818 


-2.985 


-1.233 


-1.512 


-1.880 


15 


-0.251 


-0.824 


-2.644 


-0.763 


-1.105 


-1.653 


20 


-0.012 


-0.169 




-0.395 


-0.748 


-1.486 



"our results 

''results obtained in Ref [29] 
■^results for r] = 0.2 
''results for r] = 0.4 
^results for r] = 0.6 
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